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Abstract 



It is shown that the generally accepted definition of the Poynt- 
ing vector and the energy flux vector defined by means of the en- 
ergy density of the electromagnetic field (Umov vector) lead to 
the prediction of the different results touching electromagnetic 
energy flux. The experiment shows that within the framework of 
the mentioned generally accepted definitions the Poynting vec- 
tor adequately describes the electromagnetic energy flux unlike 
the Umov vector. Therefore one can conclude that a generally 
accepted definitions of the electromagnetic energy density and 
the Poynting vector, in general, are not always compatible. 
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I. INTRODUCTION 



In the article "Motion equations of the energy in the bodies" [1] that 

appeared in the year that "Tractate" was published by Maxwell (1873), 

Umov developed the consequences from the idea of the energy localization 

in the mediums. To each volume element in the medium, the particles of 

which are in movement, an energy, constituted by the sum of the alive forces 

of the particles and elements and the potential energy, is associated. Umov 

thinks about the problem of settling down in general form "the laws of the 

transition of the energy from an element to another" , and to determinate, 

starting from general principles, the study of the movement of the energy 

in the mediums. Starting from the energy conservation law Umov deduces 

the motion equation for the energy in the mediums. If we represent the 

energy density in a given point of the medium by means of w, and through 

v x , v y and v z the energy velocity components in this point, then the energy 

density loss in that point in unit of time is determined by the relationship 

dw _ d (wv x ) | d (wv y ) | d (wv z ) 
dt dx dy dz 

"The expression (1), similar to the expression of the mass 
conservation law in hydrodynamics, is the expression of the 
elementary energy conservation law in the mediums", Umov 
writes. From this expression it can be established "the rela- 
tionship among the quantity of energy, that in unit time leaves 
toward the medium through its frontier, and the change of the 
quantity of energy in the medium" . 

This relationship is expressed with the integral expression (Umov theo- 
rem) 
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-—dxdydz + 



IS 



wv n da = 0. 



(2) 



The vector tov defines the energy flow which crosses, in the unit time, 
the perpendicular to this vector unitary surface. This is the so-called Umov 
vector. 

The case of the electromagnetic field, as particular case of the Umov 
theorem, and therefore of the Umov vector, was studied by Poynting. 

In the year 1884 J. Poynting published the article [2] that contained 
the previously mentioned Umov-Poynting theorem. In this work Poynting 
independently arrives to the same point of view developed 10 years before 
by Umov. Poynting writes: 

"If we recognize the continuity of the energy movement, that 
is to say we recognize that when the energy disappears in some 
point and appears in other, it should pass through the interme- 
diate space, then we are obliged to reach the conclusion that 
the surrounding medium contains at the least a portion of the 
energy and that it is capable to transmit the energy from one 
point to another." 

Further on Poynting, leaning on the Maxwell idea about the energy 
localization in the field, formulates in this way the main idea of his work: 

"The objective of this article is to demonstrate that there 
exists a general law for the energy transport, in agreement with 
which the energy in any point moves perpendicularly to the 
plane containing the lines of the electric and magnetic forces, 
and that the quantity of the energy passing through the unitary 
surface in this plane, for unit of time, is equal to the product 
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of the magnitudes of these two forces multiplied by the sinus of 
the angle among them and divided among 4n" . 

By this way Poynting defines the energy-flux vector for the case of the 
electromagnetic field. 

Discussing today the conception of the Poynting vector and the num- 
ber of basic difficulties associated with this concept one can sense clearly 
that neither among researchers (see, e.g., [3,4] and corresponding references 
there) nor among authors of the generally accepted text-books of classi- 
cal electrodynamics (see, e.g. [5-9]) a general agreement exists about the 
essence of the energy-flux vector related with electromagnetic fields. Actu- 
ally, the well-known authors Panofsky and Phillips state [5]: 

"Paradoxical results may be obtained if one tries to identify the 
Poynting vector with the energy flow per unit area at any particular 
point" . 

Contrarily, Feynman states [9] that exclusively the identification of the 
Poynting vector (in its generally accepted form) with the energy flow per 
unit area allows to understand the law of conservation of the angular mo- 
mentum in some special cases. Other well-known authors Landau and Lif- 
shitz state [6]: 

"Therefore the integral § Sdf must be interpreted as the flux of 
field energy across the surface bounding the given volume, so that the 
Poynting vector S is this flux density - the amount of field energy 
passing through unit area of the surface in unit time." 

Tamm [7] also identifies the Poynting vector with the energy flow per unit 
area at any particular point, however, taking into account that the definition 
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s = 

[8]: 



^(ExH) is not unique. In turn Jackson claims in his famous text-book 



"The vector S, representing energy flow, is called the Poynting 
vector. It is given by S = E x H (6.109) b... Relativistic considera- 
tions (Section 12.10) show that (6.109) is unique." 

The only way we can verify the standard formula for the energy flow 
due to the electromagnetic field is by experiment. Feynman said [9]: 

"There are, in fact, an infinite number of possibilities for w (en- 
ergy density) and S, and so far no one has thought of an experimental 
way to tell which one is right." 

In this work we theoretically rationalize that the Poynting vector (in its 
standard definition) does not always coincide with the energy flux vector 
(Umov vector) related with electromagnetic waves. The results of the ex- 
periment show that the Poynting vector is not always compatible with the 
generally accepted definition of the electromagnetic energy density. 

II. THEORETICAL MOTIVATION OF THE EXPERIMENT 

More often than not physicists implicitly suppose that the Poynting 
vector 

S = — E x B (3) 

4:71 



Jackson uses SI 



5 



and Umov (energy flux) vector 2 



U = wvn (4) 

always coincide for any electromagnetic wave spreading in vacuum in every 
point. Here n is a unit vector along the direction of propagation of the 
electromagnetic energy, v is the transferring energy velocity (in the case of 
electromagnetic waves in vacuum v = c) and w is the energy density of 
the electromagnetic wave. In actual fact, this assertion is shown at least 
for plane and spherical electromagnetic waves in vacuum (see, e.g., [6], Eq. 
47.5). Nevertheless, the assertion that S = U for waves of a more general 
kind is not proved in textbooks and monographs. 

Let us study what condition in vacuum for E and B in an electromag- 
netic wave must be satisfied when the equality S = U is valid. We have in 
CGS (Gauss' system): 

S = — E x B = — £5 sin an (5) 

47T 47T 

and 

V = wcn = —(E 2 + B 2 )n. (6) 
8% 

Equating (5) and (6) we obtain 

2EBsma = E 2 + B 2 (7) 



2 The expression of the Umov vector (U = rov) is obtained from the general 
energy conservation law = — V{u>v}) and describes the energy flux density 
of any kind of energy (not only electromagnetic energy), w is the corresponding 
energy density and v is the propagation velocity of the energy in a given point. 
Thus the Poynting and Umov vectors should always coincide. 
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or 

(E - Bf + 2EB{1 - sin a) = 0. (8) 

According to the problem definition we choose real values of E, B and a 
only, where a is the angle between E and B. Therefore the last equality 
(8) can be valid if and only if E = B and a = n/2. Thus the Theorem 
takes place: for the equality of the Poynting vector and Umov vector it is 
necessary and sufficient that E_LB and E = B. 

In the next sections we propose and perform the experiment which al- 
lows us to check the incompatibility of the conventional functional forms 
of the Poynting vector and the electromagnetic energy density when the 
electromagnetic wave field does not satisfy the conditions 

E_LB and E = B. (9) 

III. THEORETICAL PREDICTIONS 

In order to obtain theoretically the electromagnetic energy flow observed 
in the experiment described in section IV, we model the sources with the 
help of two point sources emitting spherical waves and we calculate the 
flow intensity by using the definition of the electromagnetic energy flow 
according to the Poynting vector and according to the Umov one. 

Let the source 1 be placed in (0, —I, 0) and the source 2 in (0, /, 0). The 
screen is placed in the plane z = h with —a < x < a and —b<y<b (Fig. 

1-) 
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FIGURES 



FIG. 1. Position of the point sources and the screen. 

The monochromatic spherical waves created by these sources can be 
modeled with the following expressions for the electric and magnetic fields: 

E 

E 1 = -^cos(ki.Ri-wt)e <> i, (10) 

Hi 

Bi = — - cos (ki • Ri - ut) e^i (11) 
ti\ 

E 

E 2 = cos (k 2 • R 2 - ut) e e2 , (12) 
-n-2 

B 2 = — cos (k 2 • R 2 - ut) e^ 2 , (13) 
-n-2 

where ki = ^-Ri y k 2 = ^-R2, E and B are amplitudes, k = 2ir/\ 
is wave number for both waves. By this way the energy flow has radial 
direction for each of these sources 

51 = T^Ir cos2 ( k i • R i - ut) ea, (14) 

cE B 

5 2 = , °„ 2 ° cos2 ( k 2 - R 2 - wt) e R2 , (15) 

where em, e ei , e vl and e R2 , e e2 , are the corresponding local spherical 
orts associated to the sources 1 and 2. By means of r, 9 y (p we will desig- 
nate the spherical coordinates in our coordinate system. The corresponding 
electromagnetic energy densities are defined by the expressions 



wi = E \ + D f ° cos 2 (kx ■ Ri - cut) , (16) 



w 2 = cos 2 (k 2 • R 2 - cut) . (17) 

07T./X2 

If we take into account that E = B , then for each of the two spherical 
waves, created by the sources 1 and 2, the conditions E = B and E _L B 
are fulfilled and therefore the Poynting and the Umov vectors coincide. 

51 = Ux, (18) 

5 2 = U 2 . (19) 

Let us consider now the resulting electromagnetic field by these two 
sources at the same time 

E T = Ei + E 2 , (20) 

B T = B! + B 2 . (21) 
In this case the Poynting vector and the Umov vector will have the form 
S T = ^-E T x B T = Si + S 2 + (Ex x B 2 + E 2 x B : ) (22) 

4lX 471 



and 



XJt = wren, (23) 



where 



E^ + B^ \ 
w T = = wi + w 2 + — (Ei • E 2 + Bi • B 2 ) (24) 

07T 47T 

and where n is the direction of the electromagnetic energy propagation for 
the resulting field, that is to say the unitary vector in the direction of the 
vector S r . 
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The energy flow measured experimentally is the integral over the screen 
surface (with normal k, unitary vector in the positive direction of the Z 
axis) of the temporary average of the energy flux density. For each source 
we have 



$1 = t r (Si • k) t dxdy = f P (Ui • k) t dxdy, (25) 

J —b J —a J —b J —a 



$ 2 = f [ a (S 2 • k) t dxdy = f r (U 2 • k) t dxdy, (26) 

J—bJ—a J —b J —a 

and for the resulting field according the Poynting definition 

$ P = f b f (S T -k) t dxdy (27) 

J — b J —a 

and according Umov definition 

$u= t [ a (\J T -k) t dxdy, (28) 

J — b J —a 

where the designation (. . .) t is the time average value. 
Laborious calculations show that 

chE^r r dxdy = chE B r r dxdy 

8n J J Rf 8tt J J Rl ~~ V ; 

and that the exact relationships between the resulting flow and the flows 
for separate sources (K, = ^) are 

ICp = — = 1 + (30) 
2$ 2(3 v ; 

for the Poynting definition, and for the Umov definition 

< 31 » 

where 



10 



f« cos(k(Ri - i? 2 )) 



+ R2) R2R1 , 2 



' -b J ~a U1U2R1R2 

2-(R 2 - R 1 )sm 3 9smcp 



sin^ 



- + i? 2 ) (2 sin 2 9 sin 2 y? + cos 2 9 + 2 sin 3 cos 2 y?) 



dxdy, (32) 



= 



6 /"« dxdy 



b J —a R"l 



(33) 



rb ra 



K 



: {F + GH) (J + KL) 
-bJ-aJo (A + B + C + 2D + 2E) 1 ' 2 



dtdxdy, 



(34) 



cos 4 (kRi — out) , 



(35) 



B = —j cos 4 (/ci? 2 - , 
it 2 



(36) 



\J\U2R\R2 



cos 2 (/ci?i — u;i) cos 2 (&i?2 — Cttf) 



X < 



(i?l + i? 2 ) 1 



sin 2 9 + (RxUl - R 2 Ul) l - sin 9 sin ip-{R 2 + R 1 ) ^ cos 2 



+ 



cos 2 



cos 2 v? 



(i? 2 - Ri) sin 2 sin y? + 2^ (i? 2 + #1) sin 6> cos 2 </? - ( l -j (R 2 - Ri) sin ip 



+ 



(i? 2 - Ri) sin 2 0-2- (i?! + R 2 ) cos 2 sin 9 sin y? + (i?i - i? 2 ) cos 2 



D 



R\Rl 



1 - 



cos 2 (kRi - out) cos 2 (fci2 2 - wf) (38) 
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E 



U 1 U 2 RjR 2 2 



x 



cos (kR\ — out) cos (kR 2 — out) x 

2 



— 2 cos 2 {kR\ — out) + —2 cos 2 (A;i?2 — ^) 
R\ R2 



1 + 



R1R2 n ■ 2 n , a ■ 2 n ■ 2 

— - 3 sin v + 4 sin # sin </? 



T 

r 



(39) 



1 ,2,!. D , 1 _2 



F = —3 cos — cut) + —3 cos (A;i? 2 — wt) , 



(40) 



G 



1 



\J\U2R\R2 



cos — cut) cos (A;i?2 — wt) , 



(41) 



R2R1 



H=(R 1 + R 2 ) — ^ sin 2 0-2 (i? 2 - -Ri) - sin d sin </? 



r 



(^1 + ^2) 2sin 2 #sin 2 ^ + cos 2 # + 2sin 3 #cosV , (42) 



1 ,2,l.n , 1 _2 



J = —2 cos (/ci?i — + —2 cos (A;i?2 — ^0 , 
R 1 R 2 



(43) 



if 



U1U2R1R2 



cos (/ci?i — cut) cos (kR 2 — cot) , 



(44) 



L - (1 + ^Uin 2 ^- 



!) fi + ^a_3sin 2 ^ + 4sin 2 ^sin 2 ^ + [- 
r \ r 2 J \r 



(45) 



C/i = (sin 2 9 + 2- sin 6 1 sin <p + - 



2 \ 2 



(46) 



U 2 = (sin 2 9 — 2- sin sin + - 



1 

2 \ 2 



(47) 
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R\ = r 2 + I 2 + 2lr sin 9 sin (p, 



(48) 



i?2 = r 2 + f _ 2lr gin Q gin ^ (49) 

The /Cp and /C;y values as a function of the angle a between the rays, 



were obtained numerically for the values \/l 2 + h 2 = 0.3 m, A = 632.8 nm, 
a = 3.5 mm and b = 2.5 mm used in the experiment. The graph of these 
dependencies are presented in the Fig. 2. 

FIG. 2. Theoretical curves of the behavior of the coeficient K, versus angle 
according the Poynting and Umov vectors. 

It is not difficult to show that the previous results for the relationship 
between the resulting flow and the flows for separate sources (/C = ^f^) 
are conserved if the spherical waves 1 and 2 are modeled by means of the 
equations 

Ei = — cos (ki • Ri - out) e v i, (50) 
Hi 

Bi = -— %os(ki • Ri -ut)e ei , (51) 

Hi 

E 

E 2 = — cos(k 2 ■ R 2 -cjt)e </ , 2 , (52) 
H2 

B 2 = -— cos (k 2 • R 2 - out) e e2 . (53) 
^2 



IV. DESCRIPTION OF THE EXPERIMENT 



The diagram of the experimental arrangement is shown in the Fig. 3 

FIG. 3. The experimental arrangement diagram. II and 12 are optical obtu- 
rators. ADC is the analog digital converter. 
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As coherent monochrome light source a He-Ne laser model 08181.93 from 
the company "PHYWE" with parameters: wavelength I = 632. 8nm, beam 
power P = lmW, polarization 500 : 1, beam diameter d = 0.5mm, is used. The 
laser beam goes to the beamsplitter Bl, where it is unfolded in two rays of 
the same intensity approximately. The reflected ray goes to the reference 
photodiode PIN1, model 1PP75 from the company "TESLA", that works 
in short circuit regime. The photodiode-generated current is amplified by 
an amplifier DC and by means of an analogical-digital converter ADC is 
received on the computer PC. This signal serves as a reference signal and 
gives us information about the intensity changes of the laser beam. 

The part of the laser beam, which goes through the beamsplitter Bl, 
goes to the second beamsplitter B2, where on the other hand it is unfolded in 
two rays. The second-beamsplitter-reflected ray goes to the optic obturator 
12, then reflected by the mirror M2 and falls on the lens L2. The lens L2 is a 
biconvex lens and possesses focal distance of 18mm. This lens transforms the 
cylindrical ray in a divergent beam. After going through the lens the beam 
goes to the measuring photodiode PIN2 (type 1PP75 from "TESLA") and 
it falls on its active surface under an angle a. The ray that goes through the 
beamsplitter B2, goes consecutively through the optic obturator II and the 
adjustable compensator C, then is reflected by the mirror Ml and falls on 
the lens LI. The lens LI is of the same type as the lens L2 and has the same 
function. After going through the lens LI the divergent beam, produced by 
the lens LI, falls on the measuring photodiode under the same angle. The 
angle between both beams falling on the PIN2 is equal to 2a (Fig. 3). The 
photodiode PIN2 also works in short circuit regime. Its signal is amplified 
by another amplifier DC and goes through an analogical-digital converter 
toward the computer. Its signal is proportional to the luminous flow that 
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falls on its active surface. 

The mirrors Ml and M2 are movable. The distances between each lenses 
and the measuring photodiode PIN2 are the same and they measure 30cm. 
These distances stay fixed during the experiment. The angle a is changed 
only from 14° until 86° and its value is measured within the accuracy of 
0.5°. When carrying out the experiment the result of each measurement is 
corrected with the reading of the photodiode PIN1. By this way the error 
produced by the laser instability is avoided. 

All the measurements are carried out in a dark room. To avoid the in- 
fluence of the laser instability on the experimental results, a normalization 
of the readings of the photodiode PIN2 is executed with the help of pho- 
todiode PIN1 for all the measurements. For each value of the angle a the 
experiment is executed in three stages: 

Stage 1: Both optical obturators II and 12 are closed up, and by means 
of the photodiode PIN2 the ground is measured. As the ground value was 
always below 0.5% of all other measurements, any correction is not applied 
to the experimental results. 

Stage 2: This stage has as a goal to equal and to measure the light 
energy flows that go through both optical branches in the experimental 
arrangement: branch 1 (optical obturator II, adjustable compensator C, 
mirror Ml, lens LI) and branch 2 (optic obturator 12, mirror M2, lens L2). 
The obturator II closes up and the obturator 12 opens up. The photodiode- 
PIN2-generated current is measured. Then the obturator II opens up and 
the obturator 12 closes up. By means of the adjustable compensator C the 
present current in the photodiode PIN2 is adjusted similarly to its previous 
current with an error limited to 1%. By this way the readings of the photo- 
diode PIN2 that correspond to the optical flow $! that passes through the 
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branch 1, and also correspond to the optical flow <3> 2 , that passes through 
the branch 2, are already known. 

Stage 3: Both obturators II and 12 open up and the current of the 
photodiode PIN2 that represents the total flow $ is measured. Then the 
computer PC calculates the coefficient 

$ 

K = * , (54) 

and it memorizes these values as a function of the angle a. 

The experiment was carried out for the two light flows configurations 
shown in the figs. 4-5, 

FIG. 4. Equivalent scheme of the experimental arrangement when the electric 
field vectors from the wave in the branch 1 are parallel to the electric field vectors 
from the wave in the branch 2, 2a is the angle between straights connecting each 
of two sources with the centre of the sensor 



FIG. 5. Equivalent scheme of the experimental arrangement when the mag- 
netic field vectors from the wave in the branch 1 are parallel to the magnetic field 
vectors from the wave in the branch 2. 

which correspond to the Eqs. (10)-(13) and (50)-(53). The experimental 
results are presented in the graphs (fig. 6 and fig. 7). 
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FIG. 6. The theoretical curves of the behavior of the coeficient JC versus 
angle according to the Poynting vector and according to the Umov vector and 
the experimental results for the case of the waves with parallel magnetic field 
vectors. 

FIG. 7. The theoretical curves of the behavior of the coeficient /C versus angle 
according to the Poynting vector and according to the Umov vector and the 
experimental results for the case of the waves with parallel electric field vectors. 

V. CONCLUSIONS 

Our work starts from the position that the Umov vector defines, in a 
general way, the energy flow for any type of energy and it is a consequence 
of the energy conservation law. Its expression for the particular case of 
electromagnetic waves is U = wen. On the other hand the Poynting vector 
defines the flow of the electromagnetic energy as a consequence of the energy 
conservation law in the Maxwell's theory. 

However, although both these vectors represent the same physical quan- 
tity and therefore they must coincide, we demonstrated that their equality 
is limited to the case when the fields, generating the electromagnetic energy 
flow, fulfill the conditions E = B (in CGS) and E_LB. 

Obviously these conditions are not general and therefore one can find 
situations when the electromagnetic field does not fulfill one or both these 
conditions. In order to determine which of the definitions of the energy 
density flux vector (Umov's vector or Poynting's vector), in the case when 
these conditions are not fulfilled, gives correct prediction about the energy 
flow, in the present paper we experimentally measured the flow of the resul- 
tant electromagnetic energy of two bunches of electromagnetic waves when 
aforementioned conditions are not fulfilled. This experiment shows that the 
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Umov vector does not describe appropriately the electromagnetic energy 
flow, while the Poynting vector does. 

Therefore, the experiment apparently shows that the energy flux density 
definition through the Umov vector in general is not applicable to electro- 
magnetic phenomena. However such conclusion must be necessarily incor- 
rect. Indeed, the expression for the Umov vector is obtained starting from 
the universally accepted conservation law of any type of energy. Conse- 
quently the Umov and Poynting vectors should always coincide. For this 
reason there is an apparent contradiction in the fact that the flow theo- 
retically predicted on the basis of Umov vector coincides neither with the 
experimentally measured flow nor with the corresponding flow calculated 
by means of the Poynting vector. 

The explanation for this apparent contradiction can reside in the follow- 
ing: 

Because the Umov vector is a consequence of the energy conservation 
law, the Umov vector functional dependence (U = wen) should be correct 
and therefore it is necessary to examine the elements used in its construc- 
tion, namely, w (energy density), c and n. 

The propagation velocity of the electromagnetic waves (and conse- 
quently, of the electromagnetic energy) in vacuum is unique and is equal to 
c, and therefore there is no problem. 

However, on the one hand, if the standard definition of the Poynting 
vector is correct it defines the correct direction of the electromagnetic en- 
ergy flux and, in turn, obviously, the direction of the Umov vector, then 
the standard expression of the electromagnetic density utilized to build the 
Umov vector cannot be correct. On the other hand if the standard defini- 
tion of the energy density is correct then the unit vector n (direction of the 
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Poynting vector) used for constructing the Umov vector must be incorrect 
and therefore the direction of the energy propagation is not correct. This 
would mean that the standard expression for the Poynting vector does not 
define correctly the electromagnetic energy density flux. Note, however, 
that our calculations and experiment in which we used the standard ex- 
pression of the Poynting vector do not contradict the last claim: the point 
is that the Poynting vector is defined with the accuracy of the curl of an 
arbitrary vector and for this reason, in principle, it is possible that the same 
result for the integral (27) will be obtained from another expression for the 
electromagnetic flux density. If this is the case then the direction of the 
energy flow does not have the direction of E x B. 

So for the case examined by us there is an incompatibility between the 
generally accepted definition of the electromagnetic energy density and the 
conventional definition of the energy flux density expressed by the Poynt- 
ing vector. This particular case allows us to affirm that, in general, these 
standard definitions are incompatible. 
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